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THE DEARTH OF MASSIVE, HELIUM-RICH WHITE DWARFS IN YOUNG OPEN STAR CLUSTERS 1 
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ABSTRACT 

Spectra have been obtained of 21 white dwarfs (WDs) in the direction of the young, rich open 
star cluster NGC 2099. This represents an appreciable fraction (> 30%) of the cluster's total WD 
population. The mean derived mass of the sample is 0.8 M Q - about 0.2 Mq larger than the mean 
seen among field WDs. A surprising result is that all of the NGC 2099 WDs have hydrogen-rich 
atmospheres (DAs) and none exhibit helium- rich ones (DBs), or any other spectral class. The number 
ratio in the field at the temperatures of the NGC 2099 WDs is DA/DB ~ 3.5. While the probability 
of seeing no DB WDs in NGC 2099 solely by chance is ~ 2%, if we include WDs in other open clusters 
of similar age it then becomes highly unlikely that the dearth of DB WDs in young open clusters is 
just a statistical fluctuation. We explore possible reasons for the lack of DBs in these clusters and 
conclude that the most promising scenario for the DA/DB number ratio discrepancy in young clusters 
is that hot, high-mass WDs do not develop large enough helium convection zones to allow helium to 
be brought to the surface and turn a hydrogen-rich WD into a helium-rich one. 
Subject headings: open clusters and associations: individual (NGC 2099) - white dwarfs 



1. INTRODUCTION 

It is by now well accepted that some white dwarfs 
(WDs) change their atmospheric chemical composition 
as they cool. The fact that the ratio of DA to non- 
DA stars changes as a function of effective temperature 
renders this conclusion indisputable (see e.g., Bergeron, 
Ruiz, & Leggett 1997). The distribution of WDs from the 
Sloan Digital Sky Survey (SDSS, see Figure 11 of Klein- 
man et al. 2004) clearly shows that the overall fraction 
of hot DO WDs (those showing He II lines) first peaks 
and then declines at a temperature where the DB WDs 
(those showing He I lines) begins to rise. At cooler tem- 
peratures, the DC WDs (featureless spectrum) increase 
in number in the same temperature range that the DBs 
begin to decline (the DBs being too cool to show He I 
spectral lines). However, the distribution of WDs of var- 
ious spectral types as a function of mass is much less con- 
strained. Both the DA and DB WD mass distributions 
show a primary peak at ~ 0.6 Mq (Bergeron, Leggett, 
& Ruiz 2001). Very few, if any, massive DBs have been 
found with M > 0.8M©. 

NGC 2099 is a very rich (4000 stars), young open 
star cluster in which we have found 50 WD candi- 
dates through a CFHT imaging project (Kalirai et al. 
2001a). Given the age of this cluster (650 Myrs), most of 
these WDs have temperatures of 13,000 K-18,000 K. At 
these temperatures, the number ratio of those exhibiting 
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hydrogen-rich atmospheres (DAs) to those with helium 
atmospheres (DBs) is about 3.5/1 in the field (Kleinman 
et al. 2004). In this Letter, we present spectroscopy of 
21 WDs in NGC 2099 (> 30% of the entire cluster WD 
population). The effective temperatures are such that 
we might expect a number of DB white dwarfs among 
their number, assuming the same statistics as in the field. 
Furthermore, the contemporaneous nature of these stars 
may allow us to gather some insights into the poorly un- 
derstood physical mechanisms that underlie the DA/DB 
distinction. 

2. OBSERVATIONS 

Imaging and spectroscopic observations of NGC 2099 
were obtained with the Canada- France-Hawaii (CFH), 
Gemini North, and Keck I telescopes. In our wide field 
CFHT imaging study (Kalirai et al. 2001a), we found 50 
cluster WD candidates in the central 15' of NGC 2099. 
In the companion Letter (Kalirai et al. 2005), we describe 
the current observations and data reduction. Summariz- 
ing, we obtained imaging and multi-object spectroscopy 
with both GMOS on Gemini and LRIS on Keck of three 
small fields (5. '5 x 5/5 for Gemini and 5'x 7' for Keck). 
These fields were chosen in an unbiased way to simply 
maximize the number of WD candidates for which we 
could obtain spectroscopy. The Gemini data have a high 
signal-to-noise ratio (22 1-hour exposures were obtained 
of the same field) whereas the Keck data (2 fields were 
obtained) have both higher resolution and bluer spectral 
coverage, allowing the detection of higher order Balmcr 
lines. The observations of the second field were taken 
at high airmass and so the bluest flux has been lost as 
a result of atmospheric dispersion. This, however, does 
not affect the classification of these stars as DA or DB. 

In total, we obtained spectroscopy of 24 individual WD 
candidates in the field of NGC 2099 (3 stars turned out 
not to be WDs). Therefore, despite sampling only 14% of 
the total cluster area, we include > 30% of the total WD 
population (cluster and field) given the careful position- 
ing of the fields. This is therefore the largest individual 
star cluster WD sample that has ever been spectroscop- 
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Fig. 1. — Spectra of the 21 WDs observed in the direction of NGC 2099 shown in order of decreasing brightness from top (e.g., WD1) 
to bottom (e.g., WD8), within each panel (all with 22.46 < V < 23.66). Four of the spectra are shown twice: WD19 is the same star as 
WD1, WD20 is WD2, WD22 is WD 5 and WD23 is WD6. It is probable that 4-5 of the WDs are field objects. Note that all of the WDs 
are DA spectral type. 



ically acquired. We present the 25 spectra for the 21 
WDs (4 stars are in common between the Gemini and 
Keck data) in Figure ^ Surprisingly, all stars show hy- 
drogen Balmer absorption lines and are of DA spectral 
type. A table of the derived temperatures and masses 
for most of this sample of WDs is given in Kalirai et al. 
(2005). 

3. THE SIGNIFICANCE OF OUR SAMPLE 

Given our sample of 21 DA WDs, the 3.5/1 field 
DA/DB ratio at these temperatures suggests that we 
should see ~ 4 DB WDs (considering that 4-5 of our 21 
WDs are field objects). Assuming that binomial statis- 
tics are applicable, the probability of a sample producing 
stars by chance with this expectation is ~ 1.4%. We 
can improve the statistics significantly by including spec- 
troscopy from other young open clusters. The full sample 
in a search of the literature is shown in Table 1. All of the 
turn-off stars in these clusters have masses above 2 M Q 
so that the WDs are all likely to be massive, just as in 
our sample. The total number of hot, massive WDs ob- 
served rises to 61, where 13.5 are expected to be DBs, 
and none are found. The likelihood of this happening 
purely by chance is now 0.00002%. 

Thus it appears as though massive open cluster WDs 
are always DAs, in contradiction to the expectation from 
the overall field fraction. DBs are not normally found im- 
mediately above a temperature of 30,000 K (the DB-gap) 
nor are they found below a temperature of 12,000 K at 
which point helium lines become spectroscopically invis- 



ible (Bergeron, Ruiz, & Leggett 1997). (Of course, the 
stars may still possess helium atmospheres at these cool 
temperatures). All of our stars are between these tem- 
peratures, where DBs are abundant in the field. One 
possibility is that this indicates a hitherto unappreciated 
mass dependence to the phenomenon (§EJ. Another is 
that the cluster environment may somehow lead to larger 
hydrogen envelope masses. 

4. THE ROLE OF THE CLUSTER ENVIRONMENT 

Are cluster WDs really representative of field WDs? 
The cluster environment is unlikely to give the WDs a 
common property (such as angular momentum, magnetic 
field, or rotation) as the initial conditions from the cloud 
collapse should be washed out by the huge disparity be- 
tween cluster and stellar length scales (see e.g., Menard 
& Duchene 2004). White dwarfs in a given cluster will 
all have the same metallicity; however, our sample con- 
tains 15 clusters with different metallicities ranging from 
at least Z = 0.008 to Z = 0.024. Thus, it is reason- 
able to treat the 61 WDs as individual data points - an 
assumption necessary for the statistical veracity of our 
result. 

The fact that star clusters are bound by their own grav- 
ity might lead one to imagine that they would trap gas 
lost in the advanced stages of stellar evolution and that 
this would lead to an enhanced accretion of gas (largely 
hydrogen) by the WDs in the cluster. However, this no- 
tion is problematic from the outset. Open star clusters 
have low velocity dispersions (NGC 2099 has a = 2.5 
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TABLE 1 

White Dwarfs in Young Open Clusters 
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km/s), much lower than the wind velocities from evolved 
2-4 Mq stars (~ 10-30 km/s; e.g. Barnbaum, Zucker- 
man, & Kastner 1991). As a result, the mass lost from 
the cluster is unbound and so the cluster WDs do not ac- 
crete from a gravitationally trapped reservoir of gas. The 
WDs will, nevertheless, accrete some material from the 
escaping collective mass outflow. To determine whether 
cluster WDs might have thicker hydrogen layers than 
field WDs, we consider accretion from a 'cluster wind'. 

The wind velocity (~ 20 km/s) is larger than the WD 
space velocity, so that the wind accretion rate (M = 
kit^f-p) onto a 0.85M Q WD is 



M = 5.5 x 1O~ 16 M yr' 1 



lem 



(1) 



The real unknown here is n, the local gas density due 
to the outflowing wind. Assuming the WD is located 
at a distance r from the cluster center (which is outside 
the majority of the mass-losing stars), n can be found 
from the cluster mass loss rate M c \ = 4irr 2 Vm p n, where 
V is the wind velocity. To get an expression for M c \, 
we assume a Salpeter initial mass function for the clus- 
ter from 0.1 M to 10 M . Since NGC 2099 is young 
enough that the remnant masses are small compared to 
the progenitor masses, we can also assume that all of the 
mass above the turnoff is lost. Finally, to make this a 
function of time, we assume the turnoff mass scales with 
time as T = 6 Gyr (M/M & )~ 2 - 5 . Putting this together, 
we arrive at the present day cluster mass loss rate, 



M cJ - 2.6 x 10~ 7 M Q yr 



T 



650 Myr 



-0.86 



(2) 



where we have used M c \ = 2515 M and turnoff mass M 
= 2.4 M Q (Kalirai et al. 2001a). 

Equation can now be fed back to get n and M for 
the WD as a function of cluster age. This can then be 
integrated to yield the total amount of mass accreted by 
the WD, 



The weak dependence on age is the result of the fact 
that the cluster mass loss rate is higher at early times but 
lasts for a shorter period of time. Factors of 10 in age 
correspond to about a 30% variation in accreted mass. 

The important result to note here is that equation J3J 
does not lead to more accretion than would occur in 
the field. For a young field WD, moving at 20 km/s 
through an interstellar medium of density 1 cm -3 , M ~ 
5 x 10~ 16 Af Q yr~ x . Thus, it appears unlikely that our 
observational result can be explained by an enhanced 
accretion rate for cluster WDs. Thus we turn to the 
properties of the WDs themselves for an explanation. 

5. IS THE DB PHENOMENON A FUNCTION OF WHITE 
DWARF MASS? 

If the incidence of non-DA stars is a strong function 
of WD mass, that may explain our data. In fact, the 
results of Beauchamp et al. (1996) contain a hint to this 
effect. They find that the mass distribution of DB stars 
lacks the tail to high masses that is seen in the DA mass 
distribution (and even in the smaller DBA sample). Our 
result makes this possibility much stronger statistically. 
What are the possible physical mechanisms for such a 
mass dependence? 

5.1. The Mass Dependence of Convective Mixing 

Although the chemical evolution of WDs is well known, 
the exact manner in which these transformations are 
achieved is still a matter of discussion. The appear- 
ance of DB stars at temperatures less than 30,000 K is 
thought to be due to the onset of convection in the he- 
lium layer that lies beneath the surface hydrogen layer 
(Liebert, Fontaine, & Wesemael 1987). If the hydrogen 
surface layer is thin enough (Liebert et al. considered 
very thin layers ~ lO~ 14 Af ) then it spans only a few 
pressure scale heights and convective overshooting may 
overwhelm a surface hydrogen layer. The status of this 
idea in the face of later estimates of surface hydrogen 
layer masses (Clemens 1993) is unclear. 

Although WD models do not yet exist that contain a 
detailed treatment of convective overshooting, we can at 
least map the mass dependence of the onset of convection 
in the helium layer for a variety of hydrogen surface layer 
masses, using the models of Hansen (1999). For q(H) — 
Mft/M — 10~ 7 or smaller, WDs with masses < 1 M 
will form convective zones in the subsurface helium layer 
(provided the WDs are cooler than ~ 25,000 K). How- 
ever, as the surface hydrogen mass layer becomes larger, 
there is an upper mass limit (near 0.8 M for g(H) = 
10~ 6 ) above which a convective zone will never form in 
the helium layer. Therefore, given the higher masses of 
the NGC 2099 WDs (mostly 0.7-0.9 M - see Table 1 
of Kalirai et al. 2005) compared to field WDs, the con- 
version of a DA to a DB WD through convective mixing 
may be inhibited by hydrogen layer masses that would 
still allow the transformation of 0.6 M WDs. Note, 
however, that these layers are considerably thicker than 
in the original proposal of Liebert ct al. (1987), and 
we are unaware of any quantitative model that has yet 
demonstrated the kind of mixing required in models with 
these parameters. 
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5.2. Mass Dependence of Hydrogen Removal 

It is also possible that the mechanisms that determine 
the final hydrogen surface mass are mass dependent as 
well. The most likely value of q(H) from stellar evolution 
is q(tt) ~ 1(T 4 (Iben & Tutukov 1984), a value high 
enough that it is likely to resist any attempt at convective 
mixing. The mechanisms that produce lower hydrogen 
masses are poorly understood but some are indeed likely 
to depend on mass. For example, Lawlor & MacDonald 
(2003) show that "born-again" stars will be devoid of 
hydrogen after leaving the asymptotic giant branch for 
the final time. These stars evolve through this phase as 
a result of changes in the efficiency of convective mixing, 
which may be mass dependent. Other mechanisms, such 
as the 'self-induced nova' of Iben & MacDonald (1986) 
rely on shell flashes powered by CNO burning and so may 
be more easily quenched in more massive WDs where the 
heavier elements separate out more rapidly in the higher 
gravity. The competition between nuclear burning and 
gravitational separation in massive WDs is another issue 
that needs to be investigated anew in the light of these 
results. 

5.3. Binary Evolution 

Finally, it is possible that a number of WDs in NGC 
2099 have been affected by some type of binary star evo- 
lution (interacting, merging, etc.). There is evidence to 
suggest binary evolution results in fewer DB white dwarfs 
(Liebert, Bergeron, & Holberg 2004, § 5.2), and this 
could alter the cluster DA/DB ratio (although binaries 
with non-DA components do exist, e.g., Wood & Oswalt 
1992). Hurley & Shara (2003) have recently modelled 
the WD cooling sequence of NGC 6819 (see Kalirai et al. 
2001b) using N body simulations, and they found that 
the observed morphology of the cooling sequence can be 
understood by invoking a large double degenerate and / or 
mass transfer WD binary population. These simulations 
include contributions from a number of different binary 
evolutionary scenarios, such as common envelope evolu- 
tion and mass transfer systems. Such effects could also 
play a role in the evolution of NGC 2099 WDs, which 
do not exhibit a tight sequence of points on the color- 
magnitude diagram but rather a more diffuse clump (see 
Kalirai et al. 2001a). For a 40% primordial binary frac- 
tion, Hurley & Shara (2003) estimate that 12% of all 
WDs will be double degenerates after 0.5 Gyrs (J. Hur- 



ley, 2004, private communication). However, up to 39% 
of the WDs are nonstandard single WDs (i.e., they were 
in binaries at some point and may have suffered mass 
exchange) in which evolution may have resulted in mass 
transfer and thick H envelopes. 

6. CONCLUSIONS 

Spectroscopic observations in a sample of young open 
star clusters have shown that all cluster WDs exhibit 
hydrogen-rich atmospheres. This is in contrast to WDs 
in the field, where at temperatures between 13,000 K 
and 18,000 K, ~ 25% are of DB type. The statistics are 
now good enough that the probability of this happening 
purely by chance is vanishingly small. A potential expla- 
nation is that the cluster environment might be the cul- 
prit. However, the only likely physical mechanism here - 
the accretion of H-rich gas produced by massive evolving 
stars - is found to be not higher than the rate for field ob- 
jects. The most promising scenario appears to be related 
to the fact that these cluster WDs are, on average, more 
massive than those in the field because of the young ages 
of their parent clusters. Cooling models (Hansen 1999) 
suggest that the underlying helium convection zone may 
not easily break through to the surface hydrogen layer 
in massive WDs as it does for average mass WDs. This 
scenario makes two clear testable predictions: (1) field 
DB WDs should be less massive on average than field 
DAs at temperatures below about 20,000 K, and (2) the 
field DA/DB ratio ought to be found in older clusters 
whose turnoff stars produce WDs nearer to the canoni- 
cal 0.6M n . 
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